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^ ABSTRACT 

. Aims. To probe the spectral energy distribution (SED) of the ionizing background radiation at z < 2 and to specify the 

^ ' sources contributing to the intergalactic radiation field. 

Methods. The spectrum of a bright quasar HS 1103-1-6416 (zem = 2.19) contains five successive metal-line absorption 
systems at ^abs ~ 1.1923, 1.7193, 1.8873, 1.8916, and 1.9410. The systems are optically thin and reveal multiple lines 
of different metal ions with the ionization potentials lying in the extreme ultraviolet (EUV) range (~1 Ryd to ~0.2 
■ keV). For each system, the EUV SED of the underlying ionization field is reconstructed by means of a special technique 

developed for solving the inverse problem in spectroscopy. For the Zabs ~ 1.8916 system, the analysis also involves the 
He I resonance lines of the Lyman series and the He I A504 A continuum, which are seen for the first time in any cosmic 
object except the Sun. 

' Results. From one system to another, the SED of the ionizing continuum changes significantly, indicating that the 

• ' intergalactic ionization field at 2: < 2 fiuctuates at the scale of at least ~ 0.004. This is consistent with Az < 0.01 

' estimated from He 11 and H I Lyman-a forest measurements between the redshifts 2 and 3. A radiation intensity break by 

approximately an order of magnitude at £ = 4 Ryd in SEDs restored for the Zabs = 1.1923, 1.8873, 1.8916, and 1.9410 
systems points to quasars as the main sources of the ionizing radiation. The SED variability is mostly caused by a small 
^ ' number of objects contributing at any given redshift to the ionizing background; at scales Az > 0.05, the infiuence of 

local radiation sources becomes significant. A remarkable SED restored for the 2abs = 1.7193 system, with a sharp break 
shifted to _E ~ 3.5 Ryd and a subsequent intensity decrease by ~1.5 dex, indicates a source with comparable inputs of 
both hard (active galactic nuclei, AGN) and soft (stellar) radiation components. Such a continuum can be emitted by 
(ultra) luminous infrared galaxies, many of which reveal both a strong AGN activity and intense star formation in the 
circumnuclear regions. 
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^y-^ 1. Introduction inverse problem when the shape of the underlying ionizing 



O ■ Lmes of different 10ns observed m spectra of extragalac- ^^^-^-^ proportional to the ion fractions). The simplest 

^ . tic sources indicate that the intergalactic medium (IGM) (^^^^^g) gED reconstruction involving only the ratio of the 

. . ; IS kept ionized. Fractions of individual ions are determined radiation intensities at 4 and 1 Ryd (so-called hardness of 

, by the spectral energy distribution (SED) of the ionizing ^-^^ spectrum) can be performed on the basis of the col- 

continuum; hence, a particular shape of the SED must be densities of Hell and Hi Hues detected in the Ly-a 

^ known to transform the measured line intensities into abso- ^^^^^^ ^^^^ ^^^^ measurements have been carried 

. lute element abundances. Another reason to study the SED H j jje 11 lines between redshifts z ~ 2 and 

- ■ ■ IS that It provides information about the inputs of hard (ac- _ 3 ^j^^g different lines of sight (toward HE 2347- 

tive galactic nuc ei/quasars, AGNs/QSOs) and soft (star- 4342 ^nd HS 1700+6416), and it was found that the pa- 

formmg galaxies) radiation emitters. Thus it can help to ^ameter tj =Heii/H I fluctuates from one spatial coordinate 

specify the relative contributions of different sources to the ^^^^^^^ ^^^^ ^ characteristic small scale of A, < 0.01 

integrated ionizing continuum. (Reimers et al. 1997; Zheng et al. 2004; Shull et al. 2004, 

The spectral range considered m the present paper is the 2OIO; Fechner et al. 2006, 2007). This was interpreted as a 

extreme ultraviolet (EUV) , defined as energies between the ^gg^jt ^ variable hardness of the ionizing radiation caused 

hydrogen ionization potential at 1 Ryd (13.6 eV) and -0.2 ^y diversity of the radiation sources coupled with radiation 

keV (15 Ryd), where the soft X-ray band begins. The SED transfer effects in the non-uniform IGM. The discovery of 

m the EUV range is still poorly known since direct obser- foreground quasars located transverse to the line of sight 

vations arc difficuff due to several technical and physical j^g^ ^t the positions where the Heii/Hi ratios show their 

limitations. However, the^ EUV SED can be probed indi- forest values supports the assumption that the metagalac- 
rectly by using spectral lines of different ions observed in 
the absorption spectra of distant quasars, i.e., by solving the 



continuum is restored from the measured column densities 
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tic ionizing radiation is strongly affected by local sources 
(Jakobsen et al. 2003; Worseck et al. 2007). 

A finer reconstruction of the SED requires the use of 
metal absorption lines (Cii-Civ, Siii-Siiv, Nii-Nv, Oi- 
Ovi, etc.), ionization potentials of which sample the EUV 
range of the ionizing continuum with a quite small step. 
An important fact is that at z < 2 the metal lines are the 
sole means to probe the continuum shape since with mod- 
ern facilities the He ii Ly-a forest can be studied only at 
z > 2. Different variants of the SED reconstruction proce- 
dure and its application to the analysis of specific quasar 
absorption-line systems are described in Agafonova et al. 
(2005), Levshakov et al. (2008, 2009), and Fechner (2011). 

In the present work we restore the SED of the ionizing 
background for absorption-line systems detected toward a 
bright quasar HS 1103+6416 with Zom = 2.19 (Reimers et 
al. 1995). To reconstruct the SED from metal lines, one gen- 
erally needs special absorbers: optically thin systems with 
lines of successive ions of the same element (e.g., C ii-C iii- 
Civ, Siii-Siiii-Siiv). However, they are not often found 
in quasar spectra. Thus the quasar HS 1103-1-6416 is one 
of the most favorable objects since its spectrum contains 
just several systems suitable for the SED reconstruction 
and detached from each other by ^ 0.004 — 0.17, i.e., 
by the redshift interval comparable with that in the Hen 
Ly-a forest studies. An additional advantage is that the ob- 
servational data are available in a broad wavelength range 
from 1400 A to 5700 A, which provides a large variety of 
ionic transitions and, hence, improves the accuracy of the 
results. 

The paper is structured as follows: the computational 
method used to reconstruct the continuum shape of the ion- 
ization field is described in Sect. [5] Section [3] deals with the 
observational data with attention to the consistency of the 
wavelength calibration in different parts of the spectrum. In 
Sect.m individual absorption systems are analyzed and the 
results obtained are discussed in Sect. [5j Our conclusions 
are given in Sect. [51 

2. Computational method 

The computational procedure used to reconstruct the SED 
of the background ionizing radiation from metal lines was 
considered in detail in Agafonova et al. (2005). Here we de- 
scribe it briefly to explain the basic principles of the anal- 
ysis. 

The true shape of the incident radiation can be re- 
stored from absorption systems that are optically thin in 
the Lyman continuum both in H I and He ii, i.e., from those 
with column densities 7V(H i) < 2 x 10^^ cm^^ and iV(He ii) 
< 10"'^^ cm~^ . Another requirement is that the absorbing 
gas should be in the thermal and ionization equilibrium, 
i.e., its ionization structure is determined entirely by the 
SED of the ionizing radiation and the ionization param- 
eter U = rtph/^^gas, where riph and rigas are, respectively, 
the local densities of the ionizing photons and of the gas. 
Figure [T] shows continuum shapes commonly used to ap- 
proximate the EUV range: a simple power law, an AGN 
spectrum of Mathews & Ferland (1987), and an SED of the 
intergalactic ionization field at z < 2 calculated by, e.g., 
Haardt & Madau (1996) or Fardal et al. (1998). In this fig- 
ure, tick marks indicate the ionization thresholds of ions 
usually observed in quasar absorption spectra. The column 
densities of these ions are proportional to their fractions T.^, 



which in turn are determined by the SED above the cor- 
responding ionization thresholds. Thus, having the column 
densities of ions with ionization potentials sampling some 
energy interval of the ionizing continuum, we can try to 
solve the inverse problem: namely, to reconstruct the SED 
in this interval. The procedure is iterative and consists of 
the following steps. 

Firstly, the shape of an SED should be parameterized, 
i.e., described by a set of parameters (also called factors). 
As seen from Fig. [1] in the EUV range the intensity 
can be well fitted by a simple/broken power law, i.e., by 
a set of n power-law indices at {Jn oc and breaking 

points i — I, . . . ,n. The number and definition of fac- 
tors depend on how many lines of different ions are available 
in the absorption system: it is clear that with only a cou- 
ple of lines the continuum shape cannot be reconstructed 
in great detail, i.e., the fewer the lines the fewer the fac- 
tors. With some initial guesses for the adopted factors, a 
trial SED is computed and inserted into the photoionization 
code CLOUDY (ver.07.02.01, last described by Ferland et 
al.l998). The ionization parameter U at this stage of the 
analysis is assumed constant within the absorber and is 
commonly estimated by ratio (s) of column densities of two 
or more subsequent ions of the same element (Si Ii-Si iii- 
Si IV, or C ii-C iii-C iv). If for this value of U the trial SED 
cannot reproduce the column densities of all ions observed 
in the system or some other peculiarities arise (like abun- 
dance pattern well beyond observational and theoretical 
constraints), then it is considered as inadequate and has to 
be adjusted. To evaluate the goodness of a trial continuum 
shape, we need to choose some quantitative measure ('re- 
sponse'). The choice is performed rather heuristically and 
accounts for both the information obtained with the tried 
SED and for any other a priori available information. If, for 
example, the trial SED gives a lower than the observed ratio 
of Siii/Siiv at U corresponding to the observed Cii/Civ 
ratio, then the factor values should be varied in a way to 
produce SEDs that will give a higher Siii/Siiv. In addi- 
tion, a trial SED can possibly reproduce well the observed 
column densities of silicon and carbon ions, but only at the 
expense of an extremely high relative overabundance of sil- 
icon to carbon, [Si/C] > 0.3, which is neither predicted for 
nor measured in intergalactic absorbers. In this case the ad- 
justment should go toward the continuum shapes reducing 
the [Si/C] value to the acceptable level of [Si/C] < 0.3. 
When the SED is parameterized by a single factor (simple 
power law Ji, oc J^^"), then the adjustment can easily be 
performed by a trial-and-error procedure. In case of multi- 
ple factors, the direction (in the factor space) toward the 
desired response is found by means of a special random- 
ization technique for the factor values based on the experi- 
mental design methodology (details are given in Agafonova 
et al. 2005). The adjustment of an initial SED occurs by 
stepwise movement along this direction until the SED that 
produces the required response is found. 

Now with the found SED as an input, CLOUDY is used 
to calculate the ionization curves for all ions detected in the 
absorption system. The ionization curves are then inserted 
in the Monte Carlo Inversion (MCI) code, which is a solver 
for the inverse problem in spectroscopy based on a model of 
line formation in clumpy stochastic media, i.e., in the gas 
with fluctuating density and velocity fields (full description 
is given in Levshakov et al. 2000). The MCI reconstructs 
the distributions of the gas density (and hence of the ioniza- 
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tion parameter) and the velocity along the line of sight and 
estimates the element abundances (assumed to be constant 
within the absorber) by fitting the synthetic line profiles to 
the observed ones, i.e., it works directly with the intensities 
at every point within the line profiles. Compared with the 
estimates based solely on column densities, the abundances 
obtained by the MCI are more accurate and robust, espe- 
cially in case of complex line profiles, i.e., when significant 
density and velocity gradients are expected (see Eq.(15) in 
Levshakov et al. 2000). Another advantage of the MCI is 
the treatment of blends: supposing that all lines arise in 
the same gas, we restore the underlying gas density and ve- 
locity distributions using firstly the lines with clear profiles 
and then reconstructing the profiles of blended lines on the 
basis of these distributions and ionization curves for the 
corresponding ions. The important point here is that this 
approach also allows us to check (and to correct) the even- 
tual shifts of absorption lines caused by errors in wavelength 
calibration (Levshakov et al. 2009). If with a newly found 
SED we can self-consistently describe the profiles of all ob- 
served lines, i.e., obtain ^ 1 (P^r degree of freedom) for 
individual lines along with acceptable element abundances, 
then calculations are finished. Otherwise, the next itera- 
tion of the SED adjustment is performed, eventually with 
a different response function. 

3. Observational data 

Details of the observations and data reduction are given 
in Kohler et al. (1999). Here we repeat only what is rele- 
vant for the present study. A high resolution (FWHM ~ 8 
km s""'^ ) optical spectrum of HS 1103-1-6416 was obtained 
in the wavelength coverage of 3180-5800 A with the HIRES 
spectrograph (Vogt et al. 1994) on the Keck 10 m telescope. 
Low-resolution spectra in the UV range were taken with the 
FOS and the GHRS onboard the Hubble Space Telescope 
(HST). Depending on the wavelength range, the spectral 
resolution was 0.77 A (grating G140L, AA = 1415 - 1700 
A), 1.44 A (G190H, AA = 1572-2311 A), and 2 A (G270H, 
AA = 2223 — 3277 A), which is equivalent to velocity reso- 
lutions (FWHM) of 136-163 km s^^ , 186-275 km s"! , and 
183-270 km s~^ , respectively. 

Since the data were obtained with different instruments, 
offsets in the wavelength calibration between the corre- 
sponding parts of the spectrum are possible. In our analysis 
of an absorption-line system, all transitions are fitted with 
the same model of the underlying density- velocity field and 
artificial (calibration-induced) velocity shifts between the 
absorption lines can negatively affect the final result. To re- 
veal the eventual offsets, we used several absorption systems 
with strong hydrogen and metal lines throughout the whole 
available spectral range. The velocity structure was esti- 
mated by fitting the line profiles from the high-resolution 
Keck spectra to a number of Gaussian component. The 
synthetic profiles of low-resolution lines were then calcu- 
lated by convolving the obtained velocity structure with a 
FWHM at the position of a given line and adjusting the col- 
umn density to comply with the observed line intensities. 
The calibration shifts were evaluated by cross-correlation 
of synthetic and observed profiles. 

Some examples are shown in Fig. [2] For the system at 
z = (probably a DLA), the velocity structure was es- 
timated on the basis of the doublet Can AA3439,3969 A 



in the Keck spectrum and then transferred to other ions 
assuming purely turbulent broadening (i.e., the Doppler b- 
parameters for all ions were identical). For the systems at 
z = 1.1923 and 1.8916, the lines of ions with close ionization 
potentials and close atomic weights (e.g.. Civ and Oiv, 
Sin and Sii) were selected. For the system at z = 0.7129, 
a weak Hi Lyman limit is present in the G140 part of the 
spectrum, which allows us to estimate the H i column den- 
sity as iV(Hi) (1.0 — 1.5) X 10^^ cm~^ , i.e., the system is 
optically thin. The 6-parameter for H i lines in this system 
was calculated from the ^-parameter of Mg ii doublet in the 
Keck spectrum assuming a purely thermal broadening. 

It is seen from Fig. [2] that the calibration in the G270 
part cither complies with that in the Keck spectrum (e.g., 
Mgi A2852 A at z = 0) or demonstrates a negative offset 

of 10 km s-i (Fell AA2600,2586 A at z = 0, On A834 

A at z = 1.8916). In G190, the offsets relative to the Keck 
spectrum are much larger ('-~^ 40 km s^^ ) and change only 
slowly with wavelength (Sii A764 A at z = 1.8916, Oiv 
A787 A at z = 1.1923). Contrary to this, the calibration in 
G140 is quite irregular and lines detached in the spectrum 
only by a few angstroms can be either coherent with the 
Keck spectrum (Fell A1608 A at z = 0, Hi A937 A at 
z = 0.7129) or significantly shifted (Sill A1526 A at z = 0, 
Hi A930 A at z = 0.7129). 

The absorption lines in the Keck spectrum are also not 
free from calibration errors. For example, in the system at 
z = 0, the Can A3969 A line is offset by 1.5 km s^^ relative 
to Can A3934 A, while in the system at 2 = 0.7129, the 
Mgii A2803 A line is shifted by -1.3 km s"^ relative to 
Mgii A2796 A. The calibration of lines of the same ion can 
be easily checked by a simple cross-correlation of the line 
profiles. To reveal possible calibration shifts between lines 
of different ions, several trial MCI runs are required. We 
found that the Siiv A1393 A line in the z = 1.8916 system 
was offset by 2.0 km s^-^ relative to both low (Sin, Cii) 
and high ionization lines (Civ), i.e., its shift was caused 
not by a kinematical segregation but by miscalibration of 
the wavelength scale. In general, we confirm the results of 
Griest et al. (2010) that calibration errors in Keck spectra 
are quite common and can reach a value of ^ 0.5 pixel size 
(- 2 km s-i). 

For calculations described below, all revealed calibration 
shifts were corrected, i.e., positions of lines in each absorp- 
tion system were set according to a coherent wavelength 
calibration. 



4. Analysis of individual systems 

All calculations throughout the paper are performed with 
laboratory wavelengths and oscillator strengths taken from 
Morton (2003) for A > 912 A, except for the ion Sim, 
and from Verner et al. (1994) for A < 912 A. For Sim, 
the adopted laboratory wavelength is 1206.51 A and not 
1206.50 A as Morton's (for explanation, see Sect. 2 in 
Levshakov et al. 2009). The Lyman series limit for Hei, 
Ac = 504.2593 A, is taken from the database of NIST. We 
use element abundances of the solar photosphere as given 
in Asplund et al. (2009) and Grevesse et al. (2010). The 
abundance ratio [X/Y] means \og{Nx/NY)—\og{Nx/NY)Q, 
while [X] = log(iVx) - log(A^x)0. 
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4.1. Absorption system at z = 1.9410 

The high-resolution (Keck) spectrum of HS 1103+6416 con- 
tains the H I Ly-a Une and hues of subsequent ions of both 
carbon, C ii A1334 A, C iv AA1548, 1550 A, and sihcon, Si ii 
A1260 A, Sim A1206 A, Siiv AA1393, 1402 A (Fig. |3l ob- 
served profiles are shown by black). There are absorption 
features at the expected positions of the Nv AA1238, 1242 
A lines, but their profiles are inconsistent with each other, 
probably because of blending with some forest lines. In any 
case, they can be used to set an upper limit on the nitrogen 
abundance. Clearly present in the UV part of the spectrum 
are strong lines O ill AA832, 702, 507 A, O iv AA608, 553, 554 
A, and Cm A977 A (Fig.S]). 

The column densities for metal ions from the optical 
(Keck) spectrum that can be calculated by a simple Voigt 
profile fitting are given in Table [1] The column densities 
for lines from the UV part cannot be estimated a priori 
due to low resolutions (FWHM > 200 km s~^ ) and, as 
a consequence, possible blending of many lines within one 
profile, but they can be calculated at the next stages of the 
SED reconstruction employing the gas density and velocity 
distributions derived by the MCI procedure. 

The ion fractions depend strongly on the input metal 
abundances, so that for accurate analysis we need to know 
the Hi content as well. In Fig. |21 a neighboring system, 
which is seen mainly in the C iv absorption, is clearly vis- 
ible at radial velocity v ~ 115 km s~^ , i.e., the observed 
profile of H i Ly-a is formed by a superposition of two sat- 
urated lines. Using H i lines from the low-resolution (HST) 
spectrum of HS 1103-1-6416, we deconvolve this Hi profile 
into two sub-components, with 7V(Hi) — (2.5 — 3.0) x 10^^ 
cm~^ (6 — 33.1 km s~^ ) for the z = 1.9410 system and 
N{Hi) = (0.8 - 1.1) X 10^^ cm-2 (& = 32.6 km s^^ ) for 
the neighbor (Fig. [4]). 

Firstly, simple power laws, J^, oc j^^", with the in- 
dex a ranging from 0.6 to 1.8 were tried as a continuum 
shape in the EUV range. At U determined by the mea- 
sured ratio of A^(Civ)/iV(Cii), the column densities of 
silicon ions are reproduced with power-law continua hav- 
ing indices a = 1.5 — 1.7, with the element abundances 
estimated as [C] = -0.3 - -0.4 and [Si/C] ~ 0.0 - 0.05. 
However, these continua predict small fractions of O ill, 
and one needs a rather strong overabundance of oxygen, 
[0/C] ~ 0.5, to describe the observed profiles of the Oiii 
lines. Such an overabundance is only marginally consistent 
with measurements both in metal-rich stars (metal content 
Z > 0.5^0) and in high-metallicity Hii regions that usu- 
ally show [0/C] < 0.3 (Akerman et al. 2004; Garcia- Rojas 
et al. 2004; Bcnsby & Fehzing 2006; Bresolin 2007; Esteban 
et al. 2009). That is why more complex continuum shapes 
were tried as well: namely, broken power laws defined by 
three indices and two breaking points (five factors) in the 
range 1 Ryd < < 15 Ryd (Fig. O. The boundary to 
the soft X-ray range was fixed at 0.2 keV (14.7 Ryd) and 
a = 1.8 at -E > 0.2 keV (this last setting is arbitrary since it 
does not affect the fractions of ions relevant for the present 
study) . 

At U corresponding to the measured ratio of 
iV(Civ)/Ar(Cil), the first-guess SED (dotted line in Fig.© 
underproduced the measured iV(Siiv)/Af(Siiii) (corre- 
sponding Si IV profiles are shown in the inserts in Fig. [3]). 
The latter was chosen as a response function and the search 
in the five-factor space was performed for the direction to- 



ward SEDs with a higher response. Moving along this di- 
rection in steps, we found the range of SEDs reproducing 
both the observed iV(C iv)/7V(C 11) and 7V(Si iv)/7V(Si ill). 
However, for some SEDs from this range we obtained a neg- 
ative abundance ratio [Si/C] < 0. Such an underabundance 
of silicon can occur in the gas enriched mostly by the AGB 
stars (Garnett et al. 1995; Agafonova et al. 2011), but it 
is coupled both with a high (solar to oversolar) content of 
carbon and a relative overabundance of C to oxygen, [C/0] 
> (Kingsburgh & Barlow, 1999). In the present system, 
carbon abundance for any type of the ionizing spectrum 
does not exceed half-solar value. However, to describe the 
profiles of the strong Oiii AA832, 702, 507 A lines, the rel- 
ative oxygen content should always be at least not lower 
than that of carbon, i.e., [0/C] > 0. Thus, we can exclude 
here AGB stars as a prevalent enrichment source and set 
zero as a lower limit for the relative abundance of silicon to 
carbon: [Si/C] > 0. Using this condition as a new response 
function, we performed a second iteration of the continuum 
shape adjustment. 

The range of acceptable SEDs finally obtained, i.e., 
those that allow us to describe self-consistently all observed 
profiles along with reasonable abundance ratios, is shown 
by the shadowed area in Fig. [5] The physical parameters 
derived with SEDs from this range are listed in Table[21 The 
corresponding synthetic absorption-line profiles are plotted 
by red in Figs. [3] and [H Compared to pure power laws, 
the SEDs with a break at E ^ A Ryd give a fraction of 
O III that is almost twice as large, so that now the observed 
intensities of O ill lines can be well reproduced with a mod- 
erate oxygen-to-carbon overabundance ratio, [0/C] = 0.2- 
0.3. This behavior of the Oiii fraction is easy to explain: 
the ionization potential of Oiii is 4.038 Ryd and a sharp 
decrease of the number density of the ionizing photons just 
above this energy hampers the transition from O ill to O iv, 
thus retaining the oxygen in the double-ionized stage. 

4.2. Absorption complex at z — 1.8916 

This system was described in detail by Kohler et al. (1999). 
Here we re-analyze it with the aim of restoring the contin- 
uum shape of the UV ionizing radiation. 

The available UV parts of the HS 1103+6416 spectrum 
contain clear Lyman limits for both H i and He i, making 
it possible to measure the corresponding column densitiesQ: 
N(Hi) = (2.87 + 0.03) x 10^^ cm'^ (Fig. O and iV(Hei) 
= (1.5 + 0.3) X 10^6 cm-2 (Fig. E]). We note that this is 
the first discovery of the Hei continuum in any cosmic 
object except the Sun. The spread of the 7V(Hi) values 
is determined by the uncertainty in the emission profile 
of the blend Oii+Oiii+Feiii blueward of the Hi Lyman 
edge and for iV(Hei) by the noise in the spectrum. The 
ratio A^(Hei)/A^(Hi) is sensitive to the radiation inten- 
sity ratio at 1 Ryd and 1.8 Ryd (Lyman limit for Hei) 
and in principle could be used to estimate the slope of the 
continuum spectrum in this energy range, similar to the 
way Af(Heii)/iV(Hi) is used to probe the J1/J4 ratio in 
the Hl+Heii Ly-a forest measurements. But for a given 
ionization spectrum, the value of A^(Hei)/A^(Hi) depends 
strongly on the ionization parameter, i.e., any estimation 
of this ratio can be made only by involving metal lines. 



The photoionization cross section for the He: continuum is 
taken from Samson et al. (1994). 
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The absorption lines of metal ions detected in the Keck 
spectrum are shown in Fig. |8l Essential for the SED recon- 
struction are weak but clearly present lines of Fe ii A1608 A 
(other Fe ii transitions are beyond the observed range) and 
O I A1302 A. The latter is blended with an H i line from the 
Ly-a forest, but the component that falls into the blue wing 
of H I can be deconvolved. Although the complex represents 
a superposition of several individual absorbers (marked by 
the letters A, B, C, and D in Fig. [5]), they can be separated 
if we compare the relative strengths of their C ii, C IV, Si ii, 
Si III, and Si iv lines: strong C ii, Si ii along with weak C IV, 
Si IV lines characterize the gas in the absorber A; compara- 
bly strong pairs C ii, C iv, and Si ii, Si iv are revealed in the 
absorber B; and lines of highly ionized species C iv and Si iv 
are slightly stronger than low-ionization C ii, Si ii in the ab- 
sorbers C and D. A significantly different metal content in 
the sub-systems C and D is conceivable. Additionally, there 
are superimposed absorbers seen mostly in weak C iv (in 
some cases also in Siiv) transitions (indicated by vertical 
ticks in Fig.©. 

In order to distribute between the individual absorbers 
the total column density iV(Hi) measured from the optical 
depth at the Lyman edge, we fit all available Hi lines to 
the sum of Voigt components, each of which is described 
by three parameters: the line center, the broadening b- 
parameter, and the column density. However, unlike the 
previous system at z = 1.9410, which consists of two well- 
separated absorbers with the total Af(Hi) ~ 5 x 10^^ cm~^ , 
we have here at least four close and partially overlapping 
absorbers with the total 7V(Hi) ~ 3 x 10^^ cm^^ . Coupled 
with a low resolution (FWHM ~ 200 km s^^ ) of higher 
Lyman series lines of Hi, this makes the deconvolution 
procedure numerically unstable (ill-posed), which results 
in multiple solutions (i.e., in significantly different column 
densities for the same component, depending on model, ini- 
tial guesses, exit conditions, etc.), even if the total iV(Hi) 
remains constant. 

A common approach to regularizing an ill-posed prob- 
lem is to place some a priori known constraints on the fit- 
ting parameters. In the present case such constraints can 
be derived from metal lines: the absorber B has A^(Siii) — 
2.2 X 10^3 cm-2 and iV(Siiv) = 2.1 x lO^^ cm'^ (Table 2 
in Kohler et al. 1999). The Sim A1206 A hne is saturated 
(Fig. [S]) and hence should have N{Si ill) - lO" cm'^ . The 
line is in the Ly-a forest and in principle could be blended 
with some forest feature. However, any strong Ly-a blend 
with iV(Hi) 10^^ cm~^ is excluded by the continuum 
window at the expected position of Ly-/3 A1025 A, and there 
are no candidate systems that could produce strong metal 
lines at the position of Si ill. Thus the observed saturation 
is indeed due to absorption in Sim A1206 A, although con- 
tamination with a weak, (H i) ~ 10^'^ cm~^ , hydrogen line 
(or lines) cannot be excluded (see the synthetic Si ill profile 
in Fig. [5]). The silicon ion fractions depend quite strongly on 
the gas metallicity Z, and trial calculations with CLOUDY 
have shown that for a wide range of shapes of the ionizing 
continua the required large fraction of Si ill at Si ii/Si iv 1 
can be realized only if [Z] < —0.4. Similarly, in the absorber 
D the column density of the double- ionized silicon A (Si ill) 
~ (3 - 4) X 10" cm-2 at A^(Siiv)/A^(Siii) > 3 is repro- 
duced only at low gas metalicities, [Z] < — 1.0. It follows 
from these conditions that the sub-systems B and D should 
have A(Hi) ~ 10^^ cm~^ . Compared to the absorber B, 



the absorber C reveals two times lower column densities 
of Sill and Cii, but almost similar A(Siiv) and A^(Civ) 
(Kohler et al. 1999). Assuming for C the same metallic- 
ity as in B (which looks plausible taking into account the 
closeness and similarity of line profiles), we estimate for it 
A^(H i) ^ 5 X 10^^ cm^^ . These constraints on the H i col- 
umn densities, along with setting the bounds for the com- 
ponent centers close to the position of the strongest metal 
absorption, allow us to stabilize the Hi deconvolution and 
to obtain the parameters listed in Table [31 The synthetic 
Hi profiles calculated with these parameters are shown by 
red curves in Fig. O 

The sub-system A is optically thin both in Hi, A(Hi) 
- 4 X 10^*^ cm"2 aj^j jjg ^(jje ii) - 2 x lO" cm"^ . ^j^j^ 
its great diversity of ionic transitions, this sub-system is 
suitable for the SED reconstruction. The measured column 
densities are listed in Table [TJ The red wings of several 
lines in the sub-system A are blended with correspond- 
ing components from the other sub-systems (Fig. [5]), but 
Sin and Siiv are affected only weakly and can be used 
to estimate the ionization parameter. The power-law SEDs 
with a = 0.9 — 1.2 reproduce the column densities of all 
ions with the following abundances: [C] = -0.45 - -0.55, 
[Al,Si,Fe/C] = 0.0 - 0.1, [0/C] = 0.0 - 0.3, which are consis- 
tent with standard models of chemical enrichment. The pro- 
files of both C IV A1548 A and C iv A1550 A lines come out 
slightly underestimated, but they might be blended with 
components from weak-C iv sub-systems like those seen at 
V = —120 or —90 km s~^ (tick marks in the CiV panels in 
Fig. [5]). However, these power laws also significantly under- 
produce Civ profiles in the adjacent system B, which are 
much stronger and, hence, less affected by possible blends. 
Thus pure power-law continua are probably not the best 
choice for the present absorber. 

Softer power laws {a ~ 1.5), like those considered for 
the z = 1.9410 system, reproduce all ions including Civ, 
but only at the expense of high overabundances of O and Fe 
relative to C: [C]= -0.6 - -0.7, [Si/C] = 0.15 - 0.20, [0,Fe/C] 
> 0.5. A simultaneous overabundance of both iron and a- 
elements (O, Si) to carbon at a quite high carbon content 
is neither predicted nor observed. Besides, the total column 
densities of Cii and On, A^(Cii) 3.3 x lO^-* cm^^ and 
A^(Oii) = (8.5-9.0) X lO^-* cm^2 (ggg Fig.H]), do not indi- 
cate a high oxygen overabundance. Thus, we can conclude 
that pure power-law SEDs, if adopted as the ionizing con- 
tinua, are significantly different at z = 1.8916 and 1.9410. 

The continuum shape from the previous z — 1.9410 sys- 
tem (upper boundary of the shadowed region in Fig. [5]) was 
tried as an initial guess for a broken power-law SED. With 
this SED, the Fell A1608 A and Oi A1302 A fines can also 
be reproduced with [0,Fe/C] > 0.5. In general, ion frac- 
tions of Fell and Oi are only weakly influenced by the 
continuum shape at i? > 1 Ryd and depend mostly on the 
ionization parameter U , which starts to decrease sharply at 
> ~ 4. Thus, to obtain reasonable relative abundances 
of Fe and O, we have to modify the initial continuum shape 
so that the measured ratio A(Si ii)/A^(Si iv) (along with a 
sufficient amount of C iv) will be achieved at a possible low 
U. 

The range of the broken power-law SEDs acceptable 
for the present system is shown by the shadowed area in 
Fig. [TOl with the synthetic line profiles plotted in red in 
Fig. [5] and the estimated abundances given in Table [H 
These continua are very hard, with a slope between 1 
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Ryd and 4 Ryd corresponding to a = 0.5 — 0.6 (cf. with 
a = 1.3 - 1.5 at z = 1.9410). We note that Kohler et al. 
(1999), in spite of their quite different approach, also pro- 
posed a continuum shape with a = 0.5 and an order of 
magnitude break in the intensity at 4 Ryd as a favorable 
SED for the present system. 

Such hard slopes at E > 1 Ryd are indeed observed 
in some AGNs (Scott et al. 2004; Shull et al. 2012). This 
type of objects includes a well-known quasar HE 2347-4342 
(toward which the Hen Ly-a forest was firstly resolved); 
according to Telfer et al. (2002) it should have a = 0.56. 

Now, using the H i column densities for the individual 
components (Tabled, the restored SED, and the column 
densities for metal ions from Kohler et al. (1999) to estimate 
the ionization parameters in the sub-systems A, B, C and 
D, we can calculate the total amount of He i and compare it 
with the value given by the He i Lyman limit optical depth. 
Assuming the helium abundance (in number) of 0.083, we 
obtain iV(Hei)A = (2.5 - 3.0) x lO^^ cm'^ , 7V(Hei)i3 = 
(5.0-6.0) X 10^5 cm-2,iV(Hei)c = (2.0-2.5) x lO^^ cm^^ 
and iV(Hei)£, = (2.5 - 3.0) x lO^^ cm^^ , which results in 
the total A^(Hei)tot = (1-2 - 1.5) x 10^'^ cm"^ complying 
well with the value estimated from the He i Lyman edge. 

The G140 part of the HS 1103+6416 spectrum also con- 
tains a series of He i Lyman lines. To model these lines with 
a sum of components, we need a broadening parameter for 
each component. For the component A, we obtain &Hei = 15 
km s~^ from the comparison of the Hei synthetic profile 
calculated with the density- velocity distributions estimated 
from metal lines and a single Voigt profile with the same 
A^(Hei). This means that the broadening is dominated by 
turbulent effects. Assuming the same ratio 0.9) between 
the 6-parameters of the H i and He i lines in the compo- 
nents B, C, and D, we model the He i profiles as the sum of 
four Gaussians convolved with the instrumental function. 
The synthetic He i profiles are shown by red curves in the 
right-hand panels in Fig. |9l 

4.3. Absorption system at z = 1.8873 

The system is detached from the z = 1.8916 absorber by 
only 450 km s~^ , which corresponds to ^ 2.2 Mpc if the 
offset is purely cosmological. 

The metal ions are represented by lines of C iv A1548 
A (Civ A1550 A is blended with Civ A1548 A from the 
z = 1.8916 system ) and Sim A1206 A, Siiv AA1393, 1402 
A fFig. [TT|) . There is a weak absorption at the position of 
C II A1334 A, but it is not clear whether it is indeed C ii 
or some forest feature since the spectrum of HS 1103-1-6416 
contains many weak and unidentified lines in the vicinity of 
this wavelength. In the UV part, a continuum window seen 
at the position of the O iv A553 A line and weak absorption 
at the position of Niv A765 A can be used to set upper 
limits on their column densities. The hydrogen content was 
estimated in the previous section and is given in Table [3] 
(component No. 1); the column densities for metal ions are 
Hsted in Table [T] 

Power-law SEDs with a > 1.6 reproduce the measured 
ion ratios with [C] ^ —1.0 and relative underabundance of 
silicon, [Si/C] < 0, and thus can be ruled out for the same 
reasons as in Sect. 14.11 Power laws with a = 1.6 — 1.8 give 
[C] - -1.2 and [Si/C] - - 0.2, i.e., C and Si abundances 
are in line with standard chemical enrichment. However, 



under the assumption of [0/C] ^ 0.2 — 0.3, which is also 
in line with standard enrichment, these SEDs predict quite 
a large column density of N{Oiv) — (5.0 — 6.0) x 10^"' 
cm~^ , which is only marginally consistent with the noisy 
continuum at the expected position of O IV A553 A. 

Broken power-law SEDs reconstructed for the previous 
z — 1.8916 system are definitely inconsistent with lines 
in the present absorber: at U estimated from the condi- 
tion A''(Silv) = A^(Siiii), the predicted intensity of the Cii 
A 1334 A line significantly overestimates the observed pro- 
file; besides, the measured N{Civ) requires the abundance 
ratio [Si/C] ^ —0.5. Thus, we can conclude that the SED 
aX z — 1.8873 in any case differs from that at z = 1.8916. 

Due to a small number of available lines at z — 1.8873, 
it is not possible to restrict the range of appropriate con- 
tinuum shapes as was done in Sects. 14.11 and 14.21 A large 
plurality of SEDs consistent both with the measured col- 
umn densities and with standard chemical enrichment (i.e., 
delivering [0,Si/C] > at [C] —1.0) can be proposed; 
two examples are shown in Fig. [12] by solid lines. It is seen 
that the SEDs at z = 1.8873 are considerable softer, either 
only at > 4 Ryd or in the whole EUV range, than the 
ionizing continuum at z = 1.8916. 

4.4. Absorption system at z = 1.7193 

The system is detached by 17400 km s""'^ , or ~ 90 Mpc, 
from the previous z = 1.8873 absorber. In the Keck spec- 
trum, the H I A1215 A, C iv AA1548, 1550 A, Si iii A1206 A, 
Si IV AA1393, 1402 A lines are present, with continuum win- 
dows at the expected positions of Sill A1260 A, Alii A1670 
A, and Nv A1242 A (Fig. \T^. The Cii A1334 A line is 
blended with a strong forest line, and it is unclear whether 
there is a corresponding absorption or not. The UV spec- 
trum contains a series of H i Lyman lines up to the Lyman 
limit at 2479 A (marked by the arrow in Fig. [5]). Fitting 
the Hi profiles to the three component model, we obtain 
N(Hi) = (1.5 - 1.6) X 10^^ cm-2 ^j. ^-j^g ^^^g^i hydrogen 
and 7V(Hi) = (1.4 — 1.5) x 10^^ cm~^ for the central com- 
ponent exhibiting metal lines. The column densities mea- 
sured from the Sim A1206 A, Siiv AA1393, 1402 A, and 
C IV AA1548, 1550 A profiles are given in Table [T] A highly 
unusual ratio Siiv/Civ ^ 0.5 at logU ^ — 2 (determined 
by Sim/Si IV ^ 2 and Siiv/Siii > 13) is found: in the 
systems at z = 1.9410 and z = 1.8873 with comparable 
ionization states of the absorbing gas, the Silv/CiV ratio 
is only ^ 0.06 and 0.03, respectively. The measured ratios 
Siii/Siiv, Sim/Si IV, and Siiv/C iv are neither reproduced 
by simple power-law SEDs (a = 0.5 — 1.9) nor by any of 
the broken power laws restored from the previous absorp- 
tion systems. Besides, the element abundances compliant 
with the Siiv/Civ ratio require [Si/C] > 0.6 at [C] < -2.0. 
However, at such a low [C], the overabundance of Si to C is 
well constrained by [Si/C] < 0.3 (Erni et al. 2006; Cooke 
et al. 2011). The search for direction (in five-factor space) 
toward SEDs ensuring both the measured ionic ratios and 
[Si/C] < 0.3 was performed with the same initial guess as 
for the z = 1.9410 system (Fig. \^ and with the response 
function defined as a [Si/C] value (first iteration) and a 
Siiv/Siii ratio (second iteration). The resulting range is 
shown by the shadowed area in Fig. 1121 the upper enve- 
lope gives [Si/C] ^ 0.3 and the lower one is determined 
from the condition [Si/C] ~ 0. The physical parameters 
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obtained with SEDs from the shadowed area are hsted in 
Table [2] and the synthetic spectra are shown by red curves 
in Fig. [13 

The restored SEDs differ considerably from those de- 
rived from the previous absorption systems as they are 
much softer and, most striking, have a sharp break shifted 
to lower energies. Since this shift is a decisive argument 
in the further interpretation of SEDs, it is important to 
ascertain that it is not a computational artefact. The rela- 
tive Si abundance is calculated as [Si/C] = log(Siiv/Civ) 
+ log(Tciv/Tsiiv) -log(Si/C)o, i.e., to ensure [Si/C] 
~ 0.3 the ionizing continuum should provide at U deter- 
mined by Siiii/Siiv ~ 2 at least two times lower ratio 
log(Tciv/'f Siiv) as compared to the initially tried SEDs. 
It is clear that to increase the fraction of Si iv we have to 
simultaneously push the ionization of Sim (i.e., increase 
the number of photons with E > 2.46 Ryd, the ioniza- 
tion potential of Si ill) and suppress the ionization of Si iv 
(decrease the number of photons with E > 3.32 Ryd, the 
ionization potential of Siiv). At the same time, the ioniza- 
tion of C III (the ionization potential 3.52 Ryd) should also 
be suppressed to decrease the fraction of C IV. Figure [HI 
shows photoionization cross sections for Cm, CiV, Sim, 
and Si IV as they are used in CLOUDY to calculate the 
ion fractions (T^): C ill has the highest cross section with a 
strong gradient, whereas the cross sections of Si ill and Si iv 
are smaller and more mildly sloping. This explains why the 
fractions of Cm and CiV are extremely sensitive to the 
distribution of energy in the range 3.5-4.0 Ryd, whereas 
the fractions of Sim and Siiv demonstrate a significantly 
weaker dependence. As a result of such behavior, the ratio 
U IV /Ts iiv can be lowered by two times only if the ioniz- 
ing energy is strongly suppressed starting at the ionization 
threshold of Cm, i.e., at ~ 3.5 Ryd. 

In this regard we note that the modeling of metagalactic 
ionizing continua produced by the quasar radiation trans- 
ferred through the IGM predicts at redshifts z < 2 that 
the He II breaks also shifted to lower energies (to the same 
3.5 Ryd) due to the Hubble expansion (Haard & Madau 
1996, Fardal et al. 1998). However, these breaks occur 
smeared out, and the corresponding SEDs cannot ensure 
Tfciv/Tsiiv low enough to obtain [Si/C] ^ 0.3 for the mea- 
sured Siiv and C iv column densities. For example, the ion- 
izing continuum of Haard & Madau (1996) shown in Fig.[T^ 
gives [Si/C] ~ 0.55. 

The absorber at z = 1.7193 is optically thin in the Hi 
Lyman continuum but opaque in He ii. Due to the softness 
of the ionizing continuum, the estimated iV(Heii) ranges 
from 2.5 x 10^* cm^^ (upper envelope) to 6.0 x 10^® cm~^ 
(lower envelope). This means that the incident radiation 
is attenuated within the absorber and is probably slightly 
harder at i? > 4 Ryd than the one restored from the absorp- 
tion lines. However, local processes in the absorber itself do 
not affect the energy range below 4 Ryd and thus cannot 
be responsible for the intensity suppression at 3.5 Ryd. 

4.5. Absorption system at z = 1.1923 

Moving from the z = 1.7193 system down the redshift 
scale, the next system with metal lines is found only at 
z = 1.1923, i.e., it is detached by ~ 900 Mpc (or ~ 0.3c in 
the velocity space). Some lines from this system were used 
to probe the wavelength calibration of the HS 1103-f6416 
spectrum and are shown in Fig. [5| The presence of sev- 



eral strong Lyman series lines (Hi AA1025,972 A) and a 
very shallow Lyman limit allow us to estimate the hydro- 
gen content as iV(Hi) 5 x 10^^ cm^^ . Except for the CiV 
doublet, all other expected metal lines, C ill A977 A, O iii 
AA832,702 A, Oiv A787 A, Sim A1206 A, and Siiv A1393 
A, fall in the low-resolution part of the spectrum. Their 
column densities can be estimated using the radial veloc- 
ity structure found from the analysis of the C iv lines (see 
Sect. [2]) and, hence, they are not very accurate. Besides, at 
FWHM ~ 200 km s~^ , it is always possible that some in- 
tervening absorption lines also contributed to the observed 
profile. However, with A^(Civ) = 2.2 x 10^^ cm~^ , we can 
expect the same order of magnitude column densities for 
C III, O III, and O iv, so that their possible blending with 
a weak Ly-a with A^(Hi) ^ lO^'^ cm~^ (stronger hydro- 
gen lines are excluded by continuum windows at the cor- 
responding positions of Ly-/3) will not significantly affect 
the measured column densities. Unlikely C ill, O ill, and 
OlV, lines of silicon ions with expected column densities 
of ^ 10^^ cm~^ fall in the absorption troughs and their 
column densities cannot be even roughly estimated. The 
measured column densities of 7V(0iii) ~ (4.0 — 5.0) x 10^^ 
cm-2 , iV(Ciii) ~ (1.5 - 2.0) X lO" cm^^ , and iV(Oiv) 
^ (1.0 — 1.5) X 10^^ cm~^ are sufficient to rule out both the 
power laws and soft SEDs of a type recovered at z = 1.7193 
(i.e., with a sharp break shifted to E < 3.5 Ryd) as an ion- 
izing continuum for the present system. At U determined 
by the ratio A^(Clv)/7V(Cili) > 1.0, ah these SEDs pro- 
duce a small fraction of O ill and, as a result, a large ratio 
iV(Oiv)/A^(Oiii) > 10, which significantly exceeds the ob- 
served A^(0 iv)/iV(0 III) ~ 3. The ionizing spectra with 
a break at 4 Ryd similar to those shown in Figs. [SI [TUl 
and [121 seem to be more appropriate because they provide 
iV(0iv)/7V(0iii) ~ 3-5. 



5. Discussion 

Using lines of diverse metal ions, we reconstruct the SED 
of the environmental radiation field that ionizes absorption 
systems at z = 1.9410,1.8916,1.8873,1.7193, and 1.1923. 
The SED differs significantly from system to system, so 
that we can conclude that the ionization field fiuctuates at 
the redshift scale of at least ^ 0.004. This is compa- 
rable with < 0.01 revealed from joint observations of 
He II and Hi Ly-a forest at 2.0 < z < 2.9 (Shull et al. 2004, 
2010). Figure [ni shows 77 = Ai'(Hell)/iV(Hi), the parameter 
used in the forest studies to probe fiuctuations of the ioniz- 
ing background with z. The curves r]{U) are calculated with 
CLOUDY for a range of SEDs restored for a given system 
and for the gas metallicity [Z] = —1.0 (arbitrarily chosen 
since 77- values are very weakly affected by metallicity). It 
is seen that while moving from z = 1.9410 to z = 1.8916, 
rj decreases from ~ 40 to ~ 10 (for \gU > —2 relevant for 
the majority of the forest absorbers) and then increases at 
z = 1.8874 to 30 - 50. In the vicinity of z = 1.7193, the 
77- values reach several hundreds due to a sharp decrease in 
the number of photons with energies E > 4 Ryd in the ion- 
izing continuum and, as a consequence, the impossibility to 
ionize He 11. At z = 1.1923, where the preferable SEDs re- 
turn to the range detected between z = 1.9410 and 1.8873, 
T] will be again about a few tens. The Heii/Hi forest stud- 
ies have shown a slow evolution of r/ toward smaller values 
at lower redshift, with the mean ^ ^ 70 at z = 2.9 and 
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fj AO at z ^ 2.0 (Zheng et al. 2004; Fechner et al. 2006). 
Our values are in line with these results. 



Further on, the restored continuum shapes allow us to 
identify the main sources contributing to the ionizing back- 
ground. At z = 1.9410,1.8916,1.8873, and 1.1923, these 
are definitely the AGNs. In all systems, simple power-law 
SEDs are only marginally consistent with the observed 
lines, whereas a more appropriate continuum shape cor- 
responds to a broken power law with an intensity break of 
an order of magnitude at 4 Ryd. Such an SED is usually 
attributed to absorption in the He ii Lyman continuum oc- 
curring when quasar radiation, which is assumed to have a 
simple power-law SED, passes through the IGM (Haard & 
Madau 1996, 2012; Giroux & ShuU 1997). However, study- 
ing the so-called associated systems (absorbers located in 
the physical vicinity of quasars), we found that they were 
ionized by a continuum that already had a pronounced 
break at 4 Ryd (Levshakov et al. 2008), i.e.. He ii absorption 
occurred not outside, but within the radiation-generating 
structures. A generally accepted mechanism of the majority 
of the quasar continuum is photon emission by geometri- 
cally thin but optically thick accretion disc. That optically 
thick structures are indeed involved is verified observation- 
ally by detecting the hydrogen Balmer edge in the spec- 
tra of several quasars (Kishimoto et al. 2003, 2004; Hu & 
Zhang 2012). It is important that the corresponding Lyman 
break at 1 Ryd (3.3 x 10^^ Hz) is not observed, and this 
is just what the non-LTE models of accretion discs predict 
(Hubeny et al. 2000; Blaes et al. 2001). The break is reduced 
and smeared out both by relativistic effects (Doppler shifts) 
and by summing over Hi emission and absorption edges 
from different parts of the disc. If the disc is optically thick 
in H I, it will be opaque in He ii as well, and it is clear that a 
real depth of the He ii Lyman break at 4 Ryd in the emerg- 
ing disc spectrum will be affected by the same processes 
as the H l Lyman break, primarily by a combination of the 
Hen emission and absorption. However, at energies close to 
and above 4 Ryd, the thermal accretion disc already loses 
its validity as a model for quasar continuum. At high fre- 
quencies (soft and hard X-ray bands) radiation is generated 
by other, still poorly understood, mechanisms; this means 
that calculations of the Hen emission/absorption balance 
become highly uncertain. Depending on the parameters, 
the resulting SEDs demonstrate at i? > 4 Ryd a variety of 
shapes, from a sharp flux decrease to an almost complete re- 
duction of a flux discontinuity (Hubeny et al. 2000, 2001). 
As already mentioned in Sect. [1] quasar spectra between 
~ 3 Ryd and the soft X-ray band (0.2 keV) are completely 
absorbed, so that the observational data cannot be used 
to distinguish between the model predictions. Thereby it 
seems quite plausible that the integrated spectrum of all 
radiation-generating quasar structures retains a flux break 
at 4 Ryd. This break can be further enhanced if the cen- 
tral engine is accompanied by a so-called warm absorber, 
a circumnuclear gas responsible for broad absorption fea- 
tures in the X-ray spectra and detected in ^ 50 % of AGNs 
(Piconcelli et al. 2005). With a total hydrogen column den- 
sity A^(H) ~ 10^^ — 10^'^ cm~^ , the ionization parameter 
U ^ few units, and the standard helium abundance Z — 
0.083, the predicted (by CLOUDY) column density of Hen 
is ~ 10^^ — 10^* cm~^ . Hence, it noticeably softens the radi- 
ation at E' > 4 Ryd. Thus, broken power-law SEDs restored 
for the z = 1.9410, 1.8916, 1.8873, and 1.1923 systems most 



probably represent intrinsic AGN spectra as only weakly 
attenuated, if at all, by the IGM. 

In this regard the differences between SEDs at z — 
1.9410 and z ^ 1.8916 (cf. SEDs in Figs. [5]and[ini) can be 
a consequence of different spectral shapes of the ionizing 
radiation from the contributing sources. It is well estab- 
lished that both the AGN space density and their cumula- 
tive luminosity peak around z ^ 2 (Assef et al. 2011, and 
references therein). To make some quantitative estimates, 
we employ the probability that a line of sight intersects a 
'sphere of influence' of an AGN in redshift interval Az at 
redshift z (Peebles 1993): 



which is the ratio of the physical path length 
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Here n{z) and <y{z) are, respectively, the comoving number 
density of AGNs and cross section of the sphere of influence, 
c is the speed of light, Hq is the Hubble constant, and Hm, 
f^A are dimensionless density parameters (we adopt the flat 
ACDM cosmology with Hq = 70 km s~i Mpc~\ = 0.3, 
and r^A = 0.7). 

The cross section of the sphere of influence is defined 
by a radius R at which the AGN radiation flux at 1 
Ryd, J^iR), becomes equal to the mean background meta- 
galactic flux, J^^{z). At redshift z = 1.9, the value of 
J^^ 6 X 10~^^ erg cm~^ s~^ Hz^^ sr~^ at 1 Ryd was 
measured in Tytler et al. (2004), Bohon et al. (2005), and 
Dall'Aglio et al. (2008). Taking the AGN comoving space 
density at this redshift from Wolf et al. (2003) and Groom 
et al. (2009) and choosing Az = 0.05, we calculate the in- 
tersection probability p depending on the source absolute 
magnitude Mb (blue curve in Fig. [170. The red curve m 
Fig. [TTI indicates the comoving AGN emissivity, e, at 1 Ryd. 
It is seen that bright sources with Mb < —24 (quasars) con- 
tribute the majority (~80%) of all ionizing photons, and at 
the same time their spheres of influence are the most prob- 
able candidates for intersection. Since the space number 
density of such quasars is low (at z ^ 2, n < 5 objects 
per cube of 100 Mpc side), only a small number of them 
contribute to the ionizing background at every given point 
along the line of sight. In addition, the diversity of their in- 
trinsic SEDs just causes the variations in the spectral shape 
of the resulting ionizing radiation. 

The mean free path £o for all objects with Mb < —21 
is £o = 1.5 Gpc (comoving). Thus, the probability of in- 
tersecting a sphere of influence of at least one quasar in 
the redshift interval Az - 0.05 at z - 1.9 {AX = 0.6 
Gpc) is p ~ 0.4 (it may be even higher if obscured quasars 
are taken into account). This probability is sufficient to 
suggest that at least in one of the three close systems 
{z = 1.9410, 1.8916, 1.8873) the ionizing radiation is domi- 
nated by a local source, i.e., by a quasar located transverse 
to the line of sight. 



^ Conversion of Mb into the monochromatic flux at 1 Ryd is 
performed with the spectral index a — 0.4. 
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The estimate of p can be directly tested against obser- 
vations. Searching for 'quasars near quasars', i.e., for fore- 
ground quasars focated within 30' from the hne of sight to a 
background quasar, Worseck et ah (2008) detected 7 QSOs 
per 18 lines of sight in the range 1.89 < z < 1.94, among 
them one pair of QSOs separated by Az < 0.05. We note 
here that in general such pairs are not extremely rare since 
the survey by Worseck et al. contains many other exam- 
ples at different redshifts. Thus, it is not unlikely that local 
sources affect the mean metagalactic SED even in more 
than one system between z — 1.94 and 1.89. 

The spectral shape of the ionizing radiation can also be 
altered due to absorption in the intervening systems, pri- 
marily by strong absorbers with Af(HI) > 10^^ cm^^ . The 
common assumption is that such systems are directly re- 
lated to galaxies and originate in galactic halos at different 
impact parameters. Then the mean free path between the 
halos is given by (e.g., Burbidge et al. 1977): 

4 = [fT*0*r(l + s + 2/3)]-i, (4) 

where T{x) is the gamma function. For the parameters of 
Eq. ^ the following values are taken: cr* = ttR'^ with w 
100 kpc, P « 0.35 (e.g., Matejek & Simcoe 2012), 0, ~ 
2 X 10"^ Mpc"^ (corresponding to ~ -21.5 at z 
2) and s ~ -1.25 (PoU et al. 2003; Gabasch et al. 2004). 
This provides - 8 Gpc. Given £o = 1.5 Gpc for AGNs, 
the probability of intersecting a galactic halo at an impact 
parameter < 100 kpc is several times lower as compared 
to the sphere of influence of an AGN. 

Thus, it is more likely that differences in the SEDs be- 
tween the z = 1.9410 and 1.8916 systems are indeed due to 
SEDs in QSOs but not due to the effects of radiative trans- 
fer in the IGM. We note that this is an additional argument 
supporting the intrinsic nature of the intensity break at 4 
Ryd in the outcoming AGN spectrum. 

Of course, the above consideration does not exclude that 
in some cases the shape of the ionizing background radia- 
tion can be attenuated by the intervening systems. For ex- 
ample, if we assume that the systems at z = 1.8916 and 
z = 1.8876 are illuminated by the same transverse quasar, 
then the softer SED in the latter case may indicate the 
presence of a strong absorber with iV(H I) > 5 x 10^® cm~^ 
between this system and the quasar. 

As for the SED at z = 1.7193, the footprint of a local 
source is quite certain: a deep intensity drop at E > 3.5 
Ryd points to a substantial input from a soft (galactic, i.e., 
stellar) component (see, e.g.. Fig. 1 in Giroux & Shull 1997). 
In a quasar-dominated field of metagalactic ionizing radia- 
tion, which is relatively hard at i5 > 4 Ryd, such spectral 
shapes can be sensed only in the vicinity of the source it- 
self. Possible emitters of radiation with this type of SED are 
the so-called (ultra) luminous infrared galaxies, (U)LIRGs, 
which in many cases show spectra powered by a mixture 
of AGN and starburst activity (Sanders et al. 1988). The 
(U)LIRGs that host AGNs also reveal large-scale gas out- 
flows (Alexander et al. 2010; Harrison et al. 2012). At low z 
(^ 0.05), it is even possible to resolve multiple star clusters 
with a young stellar population (Mazzarella et al. 2012) 
and superbubbles in the blowout phase, which are believed 
to form as a result of giant supernova explosions (Lipari et 
al. 2005, 2009). In addition, the z = 1.7193 absorber could 
be an outflow from a transverse (U)LIRG because the over- 
abundance of a-elements and general low metallicity of the 



gas are typical for the enrichment by massive star(s) ex- 
ploded as core-collapse supernova (SN II). 

The difference in the high-frequency range of the SED 
between pure AGN and AGN-|-starburst sources is illus- 
trated in Fig. 1181 which shows observations by Grupe et 
al. (2010) that were taken with SWIFT simultaneously in 
the optical/UV and soft X-ray bands. Blue stars indicate 
data for Mkn 335, which is a classical AGN with no star- 
forming activity, and red stars show data for Mkn 876, a 
galaxy with both strong AGN activity and intense circum- 
nuclear star formation (Shi et al. 2007; Sani et al. 2010). 
The latter SED is appreciably softer in the X-ray band. For 
comparison, SEDs restored for the z — 1.9410 and 1.7193 
systems are also shown by dotted lines. It is seen that they 
fit well in the gap between the UV and X-ray points. 

6. Summary 

The spectrum of bright quasar HS 1103+6416 (z^m = 2.19) 
available in the wavelength coverag e from 1400 A to 5700 
A contains between the redshifts z — 1.19 and 1.94 five 
successive absorption systems detached by the redshifts in- 
tervals ~ 0.004 — 0.5. These systems are optically thin 
both in H I and in He ii and reveal multiple lines of differ- 
ent metal ions, the ionization potentials of which sample 
the EUV range (1 Ryd - 0.2 keV) of the ionizing radiation 
with a fairly small step. Such properties make it possible to 
solve with sufficient accuracy the inverse problem of spec- 
troscopy: namely, to estimate the physical parameters of 
the absorbing medium simultaneously with the reconstruc- 
tion of the SED in the underlying ionizing continuum. We 
note that metal absorption systems in quasar spectra are 
the only means to probe observationally the SED of the 
ionizing radiation at redshifts z < 2. The results obtained 
are as follow. 

1. The SED changes significantly from one system to an- 
other, i.e., the intergalactic ionization field fluctuates at 
the scale of at least ~ 0.004. This is consistent with 
A^ < 0.01 estimated from Heii and Hi Lyman-a forest 
measurements between redshifts 2 and 3. 

2. In all systems, pure power law SEDs of the type J^, oc 

are not or only marginally consistent with the ob- 
served line profiles. A more preferable continuum shape 
is a broken power-law with the break at -B 4 Ryd, which 
likely represents an intrinsic SED of the outcoming AGN 
radiation. 

3. The SEDs restored for the Zabs = 1.1923, 1.8873, 1.8916, 
and 1.9410 systems point to AGNs/QSOs as the main 
source of the ionizing radiation. The SED variability is 
mostly caused by a small number of QSOs contribut- 
ing at any given redshift to the ionizing background. At 
scales Az > 0.05 the influence of local sources (fore- 
ground quasars located transverse to the line of sight) 
becomes significant. This can be tested observationally 
by a search for QSOs close to the line of sight as done 
by Worseck et al. (2008). 

4. A remarkable type of SED is restored for the Zabs = 
1.7193 system. It demonstrates a sharp break shifted to 
i? ^ 3.5 Ryd with a subsequent intensity decrease by 
^1.5 dex, suggesting comparable inputs from both hard 
(AGNs/QSOs) and soft (young stars) radiation compo- 
nents. In a quasar-dominated ionizing radiation field, 
such a soft continuum can survive only in the vicinity 
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of the radiation source itself. For the given system this 
source is probably one of the (U)LIRGs, many of which 
reveal both strong AGN activity and intense star for- 
mation in the circumnuclear regions. 
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Fig. 1. Continuum shapes used to approximate the SED in 
the EUV (1 Ryd - 0.2 keV) range: a simple power law (blue 
dotted line), an AGN continuum by Mathews & Ferland 
(1987) (blue line), an ionizing spectrum of the metagalactic 
radiation field at redshift z ^ 2 by Haardt & Madau (1996) 
(black line), and its approximation by a broken power law 
(red line). All the SEDs are normalized so that Jy{hv = 1 
Ryd) = 1. The positions of ionization thresholds of different 
ions are indicated by tick marks. 
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Fig. 2. Calibration-induced shifts of the observed profiles illustrated by spectral lines from different absorption systems. 
The velocity structure is estimated on the basis of lines from the high-resolution (Keck/HIRES) part of the HS 1103+6416 
spectrum and then applied to calculate the profiles from low-resolution (HST) parts. The synthetic profiles are shown by 
red, the vertical ticks mark positions of sub-components. indicates the velocity shift applied to the observed profile 
to align it with a reference line from the Keck spectrum. 
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Fig. 3. High-resolution hydrogen and metal absorption lines associated with the 2;abs = 1-9410 system toward 
HS 1103+6416 (black curves). The synthetic profiles calculated with the UV SEDs shown by the shadowed area in 
Fig. [S] are plotted by red curves. The zero radial velocity is fixed at z — 1.9410. The inserted figures in Si IV panels 
display the synthetic spectra (blue curves) calculated with the first-guess SED (red dotted line in Fig. [5]). 
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Fig. 4. Low-resolution lines (HST) from the Zabs = 1-9410 system. The vertical ticks in the Hi panels indicate the 
positions of individual Voigt components. The profiles of metal ions are calculated using the velocity-density distributions 
restored on the basis of high-resolution lines shown in Fig. [31 The red synthetic profiles in the O ill and O iv panels are 
calculated for two oxygen abundances [0/C] — 0.2 and 0.3 (see text for details). 



14 



I. I. Agafonova et al.: Fluctuating intergalactic ionization field at 



1.0 




Hel 


Hell 




0.5 


_ Sill 


1 

CII 


1 

Silll SilV 










' 'cui CIV 




0.0 
















NIHi oill OIV OV OVI 




-0.5 






II 1 










NIV NV 




-1.0 






\ 1 




-1.5 










-2.0 










-2.5 












0.0 0,5 1.0 1,5 



Log E (Ryd) 
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Fig. 7. He I Lyman limit at z = 1.8916 and the identified Hei Lyman series lines. The Hi Lyman limit at 2; = 0.7129 
and the Oiii A507.388 A line at z = 1.8916 are also indicated. The latter is blended with Hei Ly-^ (A507.718 A) and 
Ly-t (A507.058 A) lines. 
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Fig. 8. High-resolution lines observed in the z — 1.8916 system. Brackets indicate the sub-systems A, B, C, and D 
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Fig. 9. Low-resolution Hi and Hei lines from the z = 1.8916 system. Tick marks indicate the positions of individual 
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Fig. 10. SEDs complying with the hnes observed in the 
z = 1.8916 system: yellow cone plots simple power laws, 
Ji, oc v~", with a = 0.9-1.2; grey shadowed area indicates 
broken power laws with a break at 4 Ryd. The black dotted 
line shows an initial shape used in the SED adjustment 
procedure; here it is chosen as an upper bound of the broken 
power- law SEDs restored for the z ~ 1.9410 system (grey 
shadowed area in Fig. [S]). 
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Fig. 11. Absorption lines observed in the z = 1.8873 system. The synthetic profiles (red) are calculated with the SEDs 
shown in Fig. [T^l 
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Fig. 12. SEDs restored for the z = 1.8873 system: yellow 
cone plots simple power laws, Jy oc with a = 1.6- 

1.8; red and blue lines show two examples of broken power 
laws. The black dotted line indicates the SED restored for 
the previous z = 1.8916 system; this SED does not comply 
with the lines observed at z = 1.8873. 
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Fig. 13. Absorption lines observed in the z = 1.7193 system. The synthetic profiles (red) are calculated with SEDs from 
the shadowed area in Fig. [131 
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Fig. 14. Photoionization cross sections for the Sim, Siiv, 
C III, and C iv ions used in CLOUDY. 
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Fig. 15. Range of SEDs of the ionizing UV radiation (shad- 
owed area) reconstructed from absorption lines at Zabs = 
1.7193. The initial guess SED is plotted by the dotted line. 
For comparison, the model metagalactic SED of Haardt & 
Madau (1996) is shown by red. 
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Fig. 16. The parameter r/ = Heii/Hi as a function of the 
ionization parameter U calculated for the SEDs restored at 
z = 1.9410 (dark grey), 1.8916 (yellow), 1.8874 (light grey), 
and 1.7193 (blue). 
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Fig. 17. Blue curve: the probability for a given line of sight 
to intersect the sphere of influence of an AGN with an ab- 
solute magnitude Mb', red curve: comoving AGN emissivity 
at 1 Ryd. Both curves are calculated for z — 1.9. 



27 



I. I. Agafonova et al.: Fluctuating intergalactic ionization field at z ~ 2 




-1.0 -0..5 0.0 0.5 1.0 1.5 2.0 2.5 



Log E (Ryd) 



Fig. 18. Simultaneous UV and soft X-ray SWIFT observa- 
tions by Grupe et al. (2010) of two AGN sources Mkn 335 
and Mkn 876 (blue and red stars, respectively). Mkn 335 
is a pure AGN without significant star-formation activity, 
whereas Mkn 876 demonstrates both strong AGN activity 
and intense circumnuclear star formation. The original data 
by Grupe et al. are (i) corrected for Galactic absorption us- 
ing the X-ray cross section from Morrison & McCammon 
(1983) and iV(H) from Table 1 in Grupe et al. ; (ii) trans- 
ferred from coordinates \og{i'L^) vs. logz/ to \og{Li,) vs. 
logE; (iii) normalized by at 1 Ryd, which is found by 
linear extrapolation of the UV data to 1 Ryd with spectral 
indices reported for Mkn 335 and Mkn 876 in Shull et al. 
(2012). The solid lines in the X-ray range are the power 
laws with indices given in Table 4 in Grupe et al.. The dot- 
ted lines show the SEDs restored for the z = 1.9410 (blue) 
and 1.7193 (red) systems. 
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Table 1. Calculated column densities (in cni^^ ) for atoms and ions in systems used to reconstruct the SED of the 
underlying ionizing continuum 



xl 1 


O 11 

Cm 
Civ 


TT O T 
Ol II 

Sim Qui 
Si IV IV 


IN iV 

Nv 


A 1 TT 

Aim 


J?n TT 

re 11 






2ab. = 1-9410 








(2.5-3.0)E15 


(5.6±0.6)E12 

2.6E14° 
(2.2±0.1)E14 


(3.5±0.5)E11 (6.3-7.8)E14*' 
(8.4±0.9)E12 (1.1-1.4)1E15'' 
(1.4±0.2)E13 

zab. - 1.8916A 


<1.5E13 






(3.5-4.5)E16 


(5.2±0.2)E13 
(3.8±0.7)E12 


(8.8±0.5)E12 (9.3±0.3)E12 

9.1E13 
(2.5±0.3)E12 

2abs = 1-8873 




(7.5±0.8)E11 
3.5E11 


(2.7=b0.3)E12 


(4.8-5. 2)E15 


<3.5E12 
9.0E13° 
(7.3±0.2)E13 


(2.2±0.4)E12 (1.8-2.2)E14'' 
(2.2±0.2)E12 (3.9-4.9)E14'' 

2abs =1-7193 


<3.0E13 
<6.0E12 






(1.4-1.5)E16 


2.5E12'' 
1.0E14" 
(1.06±0.05)E13 


<4.0E11 
(9.8±0.2)E12 4.6E14'= 
(4.7±0.3)E12 2.5E14'= 

ZabB =1-1923 








~ 5E15 


(1.5-2.0)E14 
(2.2±0.2)E14 


(4.0-5. 0)E14 
(1.0-1.5)E15 









Notes. "Calculated value with C content and density-velocity distribution derived from high-resolution 
metal lines . ''Calculated value assuming [O/C] = 0.2-0.3. '^Calculated value assuming [0/C] = 0.5 



Table 2. Mean ionization parameters, C/o, and relative element abundances, [X/Y], for the systems from Table[TJ [X/Y] 
= \og{N^/NY) - log(7Vx/iVY)o. 



■Sabs 


lof 




[C/H] 


[Si/C] 


[O/C] 


[N/C] 


[Al/C] [Fc/C] 


1.9410 


-2.1 


--2.0 


-0.3 - -0.4 


0.0 - 0.1 


0.2 - 0.3 


<0.0 




1.8916A 


-3.3 


--3.2 


-0.4 - -0.5 


0.0 


0.2 




0.1 - 0.2 0.0 


1.8873 


-2.1 


--2.0 


-1.1 - -1.2 


0.0 - 0.2 


0.2 - 0.3 


<0.0 




1.7193 


-2.0 - 


- -1.75 


-2.2 - -2.3 


0.0 - 0.3 


<0.5 






1.1923 


-1.9 


--1.7 


-1.2 - -1.5 




0.2 - 0.3 







Table 3. Deconvolution parameters of the Hi profiles in the z = 1.8916 system (Fig. [5]) 



Component 
No. 


Central 
velocity V , km s~^ 


Broadening 
parameter h, km s"'^ 


Column 
density A'', cm~^ 


Comments 


1 


-446.0 


30.0 


(4.8-5.2)E15 


system z = 1.8873 


2 


-364.0 


34.0 


1.4E14 




3 


-126.0 


36.0 


9.0E13 




4 


-91.0 


14.0 


(3.0-5.0)E14 




5 


-55.4 


17.0 


(3.5-4.5)E16 


system A, Fig. |8] 


6 


9.5 


17.0 


(1.0-1. 1)E17 


system B, Fig. |5] 


7 


30.5 


15.0 


(4.5-5.0)E16 


system C, Fig. |8] 


8 


68.0 


19.0 


(8.5-9.5)E16 


system D, Fig. |5] 


9 


127.5 


25.0 


1.4E14 




10 


271.0 


29.0 


2.2E14 
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